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Q Abstract. Two-component nanoplasmas generated by strong- field ionization of 



doped helium nanodroplets are studied in a pump-probe experiment using few-cycle 

laser pulses in combination with molecular dynamics simulations. High yields of helium 
O ions and a pronounced, droplet size-dependent resonance structure in the pump-probe 

c/3 transients reveal the evolution of the dopant-induced helium nanoplasma. The pump- 

^^ probe dynamics is interpreted in terms of strong inner ionization by the pump pulse 

and resonant heating by the probe pulse which controls the final charge states detected 

via the frustration of electron-ion recombination. 
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Nanoplasmas generated by intense femtosecond laser pulses are being actively 

O investigated for gaining insights into the ultrafast djTiamics of highly excited matter on 

f—i the nanoscale, which features extraordinary characteristics. In particular the peculiar 

J> property of laser-driven nanoplasmas of emitting highly energetic particles and radiation 

K> opens up tremendous opportunities for applications as novel sources of radiation and 

for particle acceleration [H |2] . Besides, the dynamics of multi-component nanoplasmas 

turns out to crucially impact the envisaged goal of realizing single-shot ultrafast 

diffraction imaging of large natural systems in the X-ray domain [3j. 

As a result of dedicated experiments as well as simulations, the behavior of single 
component nanoplasmas in intense near infrared fs laser pulses (10^^-10^^ W/cm^) is 
fairly well- understood (see [1112] and Refs. therein). These nanoplasmas are generated 
from neutral clusters of rare-gas or metal atoms via ionization by intense near-infrared 
(IR) laser pulses. The most interesting properties of rare-gas clusters in the IR domain 
result from their high energy absorption per atom once ionized, by far exceeding the 
values achievable in atomic jets or planar solid targets [H [21 Hj. The generic picture of 
the underlying dynamics of laser-driven nanoplasmas is determined by the interrelation 
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of the nanoplasma eigenfrequency Ures and the frequency uias of the driving laser pulse. 
The dipolar eigenfrequency of the nanoplasma Ures depends on the ionic charge density 
p, which for the spherical case is hures = ^v^e^p/(3£o^e) [B E]- Resonant driving 
conditions {ures = ^las) are achieved on sub- or few-picosecond timescales when the 
plasma is diluted by ionic expansion. 

In spite of this general understanding of the plasma dynamics, experiments 
have revealed surprises in the case of nanoplasmas from two- or multi-component 
clusters: Intense IR pulse ionization of weakly doped (<1% doping) rare-gas clusters 
show enhanced electron and characteristic X-ray yields as compared their pristine 
counterparts under identical conditions [SI Ej. While in these cases the location of 
the dopant atoms could not be determined, extensive studies on the intense IR field 
ionization of metal and rare gas clusters embedded at the center of He nanodroplets 
have been carried out in the group of K.-H. Meiwes-Broer [21 [TJ El El [ID]- These studies 
highlighted the role of the embedded dopant kernel in the nanoplasma dynamics and 
ionic motion thereof. 

As compared to free metal clusters, resonance conditions were found to be reached 
at earlier delay times when measuring charge state spectra of the embedded species 
[TT1[T2]. The appearance of He"*" and He^"*" ions was related to charge transfer collisions of 
highly charged metal kernel ions with He, in accord with the observed reduced maximum 
charging of the dopants. In large droplets, caging of ionization fragments was observed, 
which induces reaggregation and the formation of snowball complexes [T2l [13] . The 
active role of the He droplet in the plasma dynamics was initially discerned by numerical 
simulations [TH [151 [13 Ej- From molecular dynamics studies is was concluded that 
the He shell undergoes rapid inner ionization and selective resonant heating following 
the avalanche-like ionization of the dopant kernel [T2l [T7] . As a consequence, well- 
separated multiple resonances as a function of pump-probe delay are predicted to occur 
in the absorption and photo electron spectra. However, no pronounced signature for 
such double resonance is seen in the pump-probe dependence of the final ion charge 
spectra of the dopant ions due to the complex interrelation of internal ionization and 
recombination dynamics. 

In the present paper we report a combined experimental and numerical study which 
unambiguously demonstrates the dynamical evolution of the host He matrix. It extends 
our recent experiments using single 10 fs pulses which have demonstrated the ignition 
of He droplets induced by only a few dopant atoms p^. Pump-probe measurements 
with 10 fs pulses performed under the same experimental conditions together with the 
corresponding numerical simulations reveal a sensitive dependence of the optimal delays 
for high charging of He ion on the size of the nanodroplet and thus the explicit role 
played by the ionized He shell in the expansion dynamics of the composite nanoplasma. 
Furthermore, signatures for multiple resonances are predicted in the He ions yields from 
the doped nanodroplets. We show that it is crucial to account for the recombination 
dynamics to capture the complete evolution of this dopant-induced nanoplasma. 
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Figure 1. Experimental arrangement: Schematic of the assembly consisting of the He 
nanodroplet source, the doping chamber, the TOF spectrometer and the beam dump. 



2. Experiment 

^He nanodroplets are an ideal tiost medium for designing well-defined nanometer-sized 
two-component clusters [19]. Owing to their weak coupling to dopant atoms and to 
their superfiuid state, He droplets can pick up other rare gas atoms in a doping cell that 
aggregate to form clusters immersed in the droplet interior [191 EO]- The experimental 
arrangement is schematically shown in Fig. [TJ A beam of He nanodroplets is produced 
by expanding pressurized ''He gas (70 - 90 bar) through a nozzle 5/im in diameter 
maintained at a temperature of 15-25 K. By varying the nozzle temperature in this 
range the mean number of He atoms per droplet is adjusted in the range 10^-10^. In a 
second vacuum chamber further downstream, the skimmed droplet beam passes through 
a doping cell consisting of a 3 cm long cylindrical cell with two collinear apertures 
(0=3 mm) to which a pressure gauge is directly attached for monitoring the cell 
pressure. By leaking into the cell a controlled amount of krypton (Kr) or xenon (Xe) 
rare gas using a dosing valve (leak rate < 10"^''mbarl/s) the mean number of dopants 
per nanodroplet K can be controlled. When taking into account the shrinkage of He 
droplets due to evaporation of He atoms induced by the pick-up and cluster aggregation 
process, K can be determined from the cell pressure and the droplet size according to a 
modified Poissonian pick-up statistics [21] ■ This model is validated using Monte-Carlo 
simulations of the pick-up process as described in Ref . [22] . 

Intense few-cycle laser pulses (~ 10 fs) at a central wavelength of 790 nm with peak 
intensities in the range 10^^ - 10^^ W/cm^ are generated by a Ti-Sapphire based mode- 
locked laser system (Femtopower, Femtolasers GmbH, Vienna). Pairs of identical pulses 
are created using a Mach-Zehnder interferometer. The delay time between the first 
"pump" and the subsequent "probe" pulse is adjusted by varying the length of one of 
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the arms of the interferometer by moving a pair of retro-reflecting mirrors mounted on 
a piezo-driven translation stage. The collinearly aligned pump and probe pulses are 
focused by a spherical mirror (focal length / = 100 mm) into the beam of doped He 
nanodroplets. Photoions are detected by a time-of-flight (TOF) spectrometer in the 
Wiley-McLaren geometry. The use of ultrashort pulses ensures that the laser couples 
purely to the electronic degrees of freedom of the (ionized) He nanodroplets whereas the 
ionic motion can be safely disregarded within the duration of the laser pulses. 

3. Theory 

For modeling the intense laser-nanodroplet interaction we employ a classical molecular 
dynamics approach [231 IE] ■ Therein ions and plasma electrons are described classically 
while atomic ionization events via tunnel and electron impact ionization are modeled 
quantum mechanically under the inclusion of local plasma field effects. Electron-ion 
recombination is taken into account for calculating ion charge state distributions. The 
key aspects relevant for the present study are sketched briefly below. 

As model systems we consider He nanodroplets {NHe=5000, 10000 and 15000) 
doped with 20 Xe atoms. The nanodroplets are initialized as spheres (fee structure) 
at the atomic density of liquid helium (p = 0.022 A~^ [IS])- After inserting the Xe2o 
core, He atoms with Xe neighbors closer than the equilibrium He-Xe distance of 4.15 A 
are removed [21]. All atoms in the doped nanodroplets are initialized in their charge 
neutral state. Electron liberation into the nanodroplet environment (inner ionization) 
via tunnel ionization (Tl) and electron impact ionization (EH) is described statistically 
via appropriate rates. The probability for TI is evaluated from the instantaneous 
Ammosov-Delone-Krainov rates [2S], employing the local electric field. A successful TI 
event results in a new plasma electron at the classical tunnel exit. Electron impact 
ionization is evaluated from the Lotz cross sections [26], taking into account local 
plasma field effects, such as depression of in-medium ionization potentials [23]. The 
resulting ions and electrons are then propagated classically in the laser field and under 
the influence of binary Coulomb interactions via 

miVi = gicEias - Vr, ^ Vij, (1) 

where rrii, qiC, and Fj are the mass, charge, and position of the i-th particle and Eias 
describes the laser electric field of two linearly polarized gaussian laser pulses. The 
pairwise Coulomb interaction V^j is described with a pseudopotential of the form 

with the elementary charge e, the interparticle distance Tjj, their charge states qi and qj, 
and a numerical smoothing parameter s. The latter regularizes the Coulomb interaction 
and offers a simple route to avoid classical recombination of electrons below the lowest 
possible quantum mechanical energy level. In our case, the smoothing parameter is 
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determined by He and has a value of s = 0.67 A. The time consuming evaluation 
of particle-particle interactions is accelerated using massive parallel computation 
techniques. 

For the identification of resonance absorption in different regions of the doped 
droplet we perform a spatially resolved analysis of the energy absorption from the laser 
field. To this end the droplet is divided into time-dependent spherical regions flxe 
and Qhc that contain the Xe core and the surrounding He droplet, respectively. The 
instantaneous power absorption in region Qk is determined from the dipole velocity by 

P^.^t)= J2 eg.r^-Eias(t), (3) 

leading to the accumulated energy absorption 

W^atW= r P'''{t')dt'. (4) 

For a meaningful comparison of experimental and theoretical results, final ion 
charge spectra have to be determined from the simulation. Recent studies have shown 
that recombination has a crucial impact on the charge distributions and thus must be 
taken into account for predicting realistic charge spectra [16], |l7l [23] . The two relevant 
mechanisms for recombination of quasi-free electrons with atomic ions after the laser 
excitation are radiative recombination and three-body recombination (TBR). As has 
been estimated previously ^3], radiative recombination can be neglected due to low 
rates [27]. TBR, i.e. electron capture after the collision of two quasi-free electrons 
in the vicinity of an ion, proceeds mainly to high Rydberg states of the ion and can 
therefore be treated classically. Hence, TBR is automatically included in reasonable 
approximation within the classical MD propagation. Besides, the dominant contribution 
of TBR, even higher order collisional recombination processes (four-body, five-body, etc.) 
are accounted for because of the fully microscopic description of classical particle-particle 
correlations. The key advantage of the direct microscopic treatment of recombination 
is that no approximations such as quasi charge-neutrality or a thermal electron velocity 
distribution need to be used. In addition, the local field effects due to screening and 
potentials of neighboring ions are included. To approximate the resulting charge spectra 
from the MD simulation, electrons are treated as recombined when bound to a specific 
ion after Ips of propagation subsequent to the probe pulse. For details see [T7] . 

4. Results 

4.I. Pump-probe dynamics 

Fig. [2] shows a typical example of the measured He^ and He^"*" ion yields as a function 
of pump-probe delay. Both ion signals exhibit a pronounced pump-probe dynamics on 
the sub-ps time scale with maximum yields for Topt ~ 0.45 ps delay. Note that the He^"*" 
signal shows higher pump-probe contrast and reaches a maximum at a slightly smaller 
value of the optical delay {Aropt ~ 20fs). 
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Figure 2. Comparison of He+ and He^"*" ion yields as a function of pump-probe delay 
when Xei5@Hei5ooo nanodroplets are exposed to two identical pulses (~ 10 fs) of peak 
intensity 7 x 10^'' W/cm^. The He^"*" yield is maximized at a delay time of 457 (±15) fs 
and the He''" yield at 476 (±15) fs. Lines are to guide the eye. 



In a simplified picture, the existence of an optimal delay can be interpreted as a 
signature for resonantly enhanced charging pUl [T21 HTj. After the pump pulse ionizes 
the doped droplet and creates an overdense nanoplasma, the probe pulse excited the 
system resonantly after the appropriate degree of expansion. We recently reported a 
high degree of ionization in the He shell induced by a few dopant atoms at the center of 
the droplet [18]. At the initial near-solid atomic densities prevalent within the droplet 
during the pump pulse, the dipolar eigenfrequency {hwres ~ 3.2 eV) by far exceeds the 
frequency hwias ~ 1.6 eV of the laser pulses at the wave length of 790 nm. Hence, 
ionic expansion via Coulombic and hydrodynamic forces is required to reach resonant 
conditions. As a consequence, the absorption of energy from the probe laser pulse by 
the ionized droplet rises sharply, which is likely to enhance electron emission as well as 
electron-ion recombination taking place in the later stage of the expansion. A detailed 
analysis of the full pump-probe evolution of the nanoplasma is presented in the following 
section. 

The delay dependence of He ion yields are investigated for various experimental 
parameters including He droplet sizes for two different dopant species - Xe and Kr. 
The He"*" and He^+ resonance curves peaked around 500 fs turns out to be a robust 
feature with respect to variations of the laser pulse intensities as well as of the species 
and number of dopant atoms. The most pronounced peak shifting and broadening is 
observed when the size of the He droplets is varied. Fig. [3] presents the delay dependence 
of He^"*" ion yields for the cases of Xe and Kr doping. The mean number of dopants per 
droplet (~ 15) is nearly the same in both cases. Clearly, the optimal delay Topt increases 
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Figure 3. Delay dependence of He^+ ion yields for various droplet sizes (mean number 
of He atoms per droplet is indicated in the legend) for the case of Xe doping. The 
peak intensity of the pump and probe pulses is 7 x 10^^ W/cm^. The mean number of 
doped atoms in the droplets is 15 ±3 in the case of Xe. The data points are fitted with 
a fifth-order polynomial from which the optimal delay values (Topt) are extracted. 
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Figure 4. (a) Optimal delay times for Xe, Kr and Ar doped He nanodroplets of various 
sizes. The Xe and Kr cases are the same as shown in figure [3] The corresponding 
variation of the widths (FWHM) of ion yield curves is shown in panel (b). The mean 
number of doped Xe and Kr atoms is 15 ± 3 and 14 ± 2.8, respectively. The peak 
intensity of the pump and probe pulses is 7 x 10^^ W/cm^. 



with increasing droplet sizes. This is similar to the delay dependence measurements of 
optical absorption by pure Xe clusters of different sizes observed by Zweiback et al. that 
were rationalized mainly by geometrical effects [28]. Figs. |4[a) and (b) illustrate the 
measured size dependences of the optimal delay and peak width (FWHM) as obtained 
from fitting the pump-probe data (Fig. [3]) with a fifth-order polynomial. Both quantities 
are found to vary roughly linearly with the average number of He atoms per droplet in 
the considered size range. 
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4-2. Simulation results 

In order to extract the mechanisms underlying the experimentally observed pump-probe 
dynamics we have performed MD simulations with doped He nanodroplets Xe2oHeAr 
exposed to two 10 fs pulses at the intensity J = 7 x 10^^ W/cm . Note that this intensity 
is well below the threshold for barrier suppression ionization of He but sufficient for 
tunnel ionization of Xe [29] • 

The simulations show that, irrespective of the nanodroplet size, charging begins in 
the pump pulse with tunnel ionization of the dopant atoms residing at the droplet center. 
Laser heating of the first released electrons induces an impact ionization avalanche that 
quickly increases the charge states of the dopant core and the closest He shells around 
it. Subsequently, triggered by the plasma produced in the core region, the rest of the He 
nanodroplet becomes ionized from inside-out and expands. Fig. [5] shows the calculated 
delay dependence of key observables for three different droplet sizes, N = 5000, 10000, 
15000, from left to right. 

Focusing on the smallest nanodroplet first, the following general conclusions can 
be extracted (see left panels in Fig. [s]). The average inner charge state of He, i.e. the 
average number of electrons removed from He atoms, is significantly enhanced for short 
pulse delays. This refiects the density dependence of the inner ionization rates, as both 
the probability for electron-ion collisions as well as the infiuence of local plasma fields 
decrease with decreasing density. Note that no direct signature of resonant heating can 
be found in the pump-probe traces of inner ionization. In contrast to that, pronounced 
maxima occur in the energy absorption of the He nanodroplet (around 175 fs) and the 
absorption in the Xe region (around 450 fs), see Fig. [sFa). The number of continuum 
electrons (outer ionization) exhibits a similar dynamics as the total energy absorption, 
purple curve in Fig. [std). The maxima in absorption and outer ionization clearly 
indicates resonant heating of the nanoplasma within the droplet. 

For a realistic comparison of theory and experiment, charge spectra have been 
determined by taking recombination processes into account. The orange line in Fig. |5](d) 
indicates the calculated final average charge states of the He atoms after applying 
our simplified recombination scheme. The effect of recombination is reflected by the 
difference between average inner ionization and flnal charge states. In contrast to the 
inner charge state, the final charge state shows a strong delay dependence over the whole 
investigated range including a pronounced maximum for a delay of 150 fs and a shoulder 
near At =450 fs. Obviously, the dynamics in the final charge states is a result of the 
interplay of inner ionization and recombination. 

Because of the high temperature dependence of three-body-recombination, 
recombination is substantially suppressed in the case of enhanced heating of the 
nanoplasma. Comparison of inner and final charge states with the absorption in Fig. [sta) 
shows that recombination is ineffective for delays with high absorption. This trend is 
most pronounced for resonant heating of the He shell. Furthermore, also the resonant 
heating of the Xe core is sufficient to notably suppress recombination and leads to 
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Figure 5. Calculated pump probe dynamics of Xe2oHejv after exciation with two 
10 fs pulses (intensity: / = 7 x 10^'*W/cm^) for three different droplet sizes, N=500G, 
10000, 15000 from left to right. The upper panels (a-c) show the energy absorption in 
the He shell and in the Xe region (as indicated) as a function of pulse delay. In the 
bottom panels (d-f ) average charge states for inner and outer ionization as well as final 
average charge states including recombination are given for different time delays. 



the shoulder in the final charge around 450 ps delay. The peak-structure in the final 
charge states can be traced back to the combined action of inner ionization and resonant 
heating, in agreement with previous results [T7] . 

A comparison of the results for the different droplet sizes reveals the same trends 
and shows that the above picture is generic for the pump-probe excitation of the doped 
nanodroplets and is in reasonable agreement with the experimental findings. For all 
sizes of the He matrix, a pronounced peak in the final charge state is predicted. The 
peak values appear for delay times similar to the experiment and are shifted towards 
longer delays with increasing size of the He droplet. The increase of the optimal delays 
can be traced back to the fact that the larger matrix requires more time to expand to 
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resonant conditions. 

5. Conclusions 

We have studied the evolution of dopant-induced He nanoplasmas with a combination 
of experiment and classical MD simulations. High yields of He^ and He^"*" ions are 
experimentally observed when doping He nanodroplets with a few Kr or Xe atoms 
upon irradiation with pairs of few-cycle IR laser pulses. The He ion signals exhibit 
a pronounced resonance feature as a function of the pump-probe delay time, which 
most sensitively depends on the average size of the He droplets. A linear increase in 
the optimal pump-probe delay as a function of the size of the nanodroplets clearly 
elucidates the dynamical role played by ionized He shells in the resonance of the 
composite nanoplasma. Our detailed numerical studies lead to the development of 
the following dynamical picture of the two-component nanoplasma evolution: (i) 
The whole interaction dynamics of the doped He nanodroplets is launched by tunnel 
ionization of the dopant atoms. Heating of the first released electrons induces an impact 
ionization avalanche that charges up the whole droplet from inside to outside, (ii) The 
observed delay dependence of the charge states can be clearly attributed to the varying 
recombination efficiency, which is minimal for resonant heating conditions, (iii) The 
optimal delay for high ion yields depends on the droplet size and the laser intensity 
(see Ref. [301). Therefore the experimental data are the result of a convolution of 
the droplet size distribution and the intensity profile of the laser focus. This leads to 
broader structures and masks details such as the numerically observed sideband due to 
the nanoplasma resonance at the core. Nevertheless, the salient features, i. e. the shift 
and the broadening of the observed resonance structures as a function of droplet size, 
are reproduced in the simulations. 

The present studies motivate further experimental and theoretical investigations of 
the dynamics in two-component nanoclusters. While the present study brings to light 
the active role of the He shells in the nanoplasma expansion dynamics, it also opens up 
interesting questions for further investigation such as the effect of the location of the 
dopant atoms within the droplet (center versus surface) on the nanoplasma ionization 
and expansion dynamics. 
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